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^ ■ ABSTRACT 

<: 

Q , We investigate the effect of a magnetic field on the global oscillation modes 

[ of a rotating fluid star in the magnetohydrodynamic approximation. We present 

^ ! general equations for the modification of any type of fluid mode due to a general 

^ ! magnetic field which is not aligned with the star's spin axis. In the case of any 

! internal dipole magnetic field we derive the equations for the frequency correc- 

tions to the r-modes. We solve for the frequency correction explicitly for the case 
when the internal dipole field is force-free, including the uniform density case. 
In the weak-field limit, the spatial form of the r-mode velocity perturbation is 
^ I unchanged, but the magnitude of the frequency in the rotating frame increases. 
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to '. neutron 
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1. Introduction 

The observation of oscillations of the Sun and other stars has led to great advances in 
our understanding of the internal structure of stars. While the oscillations of non-rotating, 
non-magnetic stars are well understood, less is known about the oscillations of stars which 
are rotating and possess a magnetic field (Cox 1980; Unno et al. 1989). The main classes 
of stars for which magnetic fields and rotation both play important roles are the rapidly 
oscillating Ap (RoAp) stars, the magnetic white dwarf stars and neutron stars. The Ap 
stars are hydrogen burning stars with strong magnetic fields. Rapid oscillations (relative 
to the spin frequency) have been observed in the subclass of stars known as RoAp stars 
(Kurtz 1990). These oscillations have been identified as high overtone p-modes. In addition, 
evidence for line splitting by rotation has been found. These stars are modeled with a 
magnetic field tilted with respect to the spin axis of the star. 
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Numerous white dwarf stars have been observed to oscillate (Gautschy & Saio 1996). 
The oscillations of the white dwarfs have been identified as g-modes. So far no magnetic field 
larger than ~ kG (Schmidt & Grauer 1997) has been observed for any of the pulsating white 
dwarfs. On the other hand, a class of white dwarfs with magnetic fields up to lO^G (Wickra- 
masinghe & Ferrario 2000) known as magnetic white dwarfs (MWD) has shown no evidence 
of pulsations. The lack of pulsations in the MWDs may be due to magnetic suppression of 
pulsations or to pulsation amplitudes which are not large enough to be detectable. 

The observation of quasiperiodic pulsar microstructurc (Boriakoff 1976; Cordcs ct al 
1990) has led to the suggestion (Van Horn 1980) that neutron star oscillations had been 
observed. While the evidence for electromagnetic detection of neutron star oscillation modes 
is inconclusive, it may be possible that neutron star oscillations will be detectable by ad- 
vanced gravitational wave detectors. It has also been proposed that the soft gamma repeaters 
(SGR) are high magnetic field neutron stars with fields as high as lO^^G (Thompson & Dun- 
can 1995). This hypothesis is supported by observations of objects such as SGR 1900+14 
which appears to have a magnetic dipole field strength in the range ~ 2 — 8 x lO^^G (Kou- 
veliotou 1998). Torsion modes in the crusts of these neutron stars may have been detected 
in SGR outbursts (Duncan 1998). 

The relative importance of magnetic fields and rotation in different types of stellar 
systems varies by many orders of magnitude, as can be seen in Table 1, where the ratio 
of magnetic field energy, Ai = to rotational kinetic energy, T = Af7?^(27r/P)^/5, 

is given for typical stars. In Table 1 we give values for the average magnetic field, and 
the magnetic field energy is calculated as though the star's magnetic field is uniform. For 
millisecond pulsars and typical white dwarf stars, rotation is much more important than 
magnetic fields at global scales. At the other end of the scale, if the soft gamma repeaters 
can be modeled as high magnetic field neutron stars, magnetic fields are dominant over 
rotation. The magnetic white dwarf stars and some high magnetic field neutron stars have 
both magnetic field and rotation energies of similar sizes. Clearly any mathematical method 
which treats the ratio of magnetic field energy to rotational energy as a small perturbation 
can not describe the phenomenology of all types of stars listed in Table 1. 

In this paper we introduce a formalism for computing global oscillation modes of mag- 
netic, rotating Newtonian fluid stars. Our method makes use of a slow-rotation expansion, 
but all other aspects are nonperturbative, in that the modes of stars with any value of the 
ratio of M./T can be computed. The method of computation is particularly useful for the 
computation of magnetically modified r-modes and inertial modes which vanish in the limit 
of zero rotation. For illustrative purposes the magnetic corrections to the r-modes of a con- 
stant density star with a homogeneous magnetic field inclined at arbitrary angle to the star's 
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spin axis will be presented. In the example presented, the computed modes are only valid in 
the limit of small M./T. The full treatment of the modes of stars with strong magnetic fields 
is more involved than our example and will be presented elsewhere. It is straight-forward 
to extend these calculations to stars with realistic equations of state and more complicated 
magnetic fields. 

There is already an extensive literature on the pulsations of magnetic stars. Radial 
oscillations of nonrotating magnetic stars were first investigated by Chandrasekhar & Fermi 
(1953). Ledoux & Simon (1957) computed the magnetic corrections to the nonradial f- 
modcs of a nonrotating star. In these early calculations, the overall effect of the magnetic 
field was assumed to be weak, in the sense that the energy stored in the magnetic field is 
small compared to the star's gravitational binding energy. While this is generally true for 
all stars of interest, the magnetic pressure can dominate over gas pressure near the star's 
surface. This suggests that a perturbative approach will break down near the surface of a 
star. Carroll et al. (1986) introduced a cylindrical model of the region near the magnetic 
pole at the surface of a neutron star and computed the spectrum of nonradial oscillations of 
the model problem. A magnetic boundary layer approach was adopted by Dziembowski & 
Goode (1996) and Bigot et al. (2000) as an alternative method for finding the magnetically 
modified modes near the stellar surface. 

Computations of nonradial oscillations of rotating magnetic stars are more complicated 
than the calculation of modes of magnetic nonrotating stars. One approach is perturbative: 
the star's rotation and magnetic field can be taken to be small corrections to the modes 
of nonrotating, nonmagnetic stars. This perturbation approach has been applied to the 
calculation of magnetically modified pressure modes of rotating stars (Shibahashi & Takata 
1993; Takata & Shibahashi 1994, 1995). The work on modes of rotating magnetic stars 
has so far focused on the pressure modes. However, in the absence of a magnetic field 
a rotating star has another family of modes, the inertial modes, which are driven by the 
Coriolis force. Malkus (1967) computed the magnetic generalization of these modes for the 
case of a star whose magnetic field is generated by an electric current which is parallel 
to the fiuid's angular velocity. In this special case the equations describing normal modes 
can be separated in oblate spheroidal coordinates. It is of interest, however, to understand 
the magnetic generalizations of the inertial modes for magnetic fields which are tilted with 
respect to the star's rotation axis, as is the situation common in astrophysics. If excited 
in the Ap stars, these modes may create beat frequencies with the high order p-modes 
which may be observable. (The inertial modes have a very small radial component to their 
oscillation and are unlikely to provide an observable luminosity variation on their own.) 
Another motivation for computing the magnetically modified inertial modes is due to the 
prediction that it may be possible for a subset of the inertial modes, the r-modes, to be driven 
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unstable by gravitational radiation (Andersson 1998; Friedman & Morsink 1998; Lindblom 
et al. 1998) in young neutron stars. In a series of papers Rezzolla et al. (2000, 2001a,b) 
have examined the consequences of a large amplitude r-mode and shown that it would have 
the effect of amplifying the star's magnetic field and hindering the further growth of the 
mode. This effect could play an important role in the early evolution of a neutron star's 
spin and magnetic field. Further work by Ho & Lai (2000) provided further insight on the 
effect of a growing r-mode on the magnetic field and calculated the form of magnetically 
modified r-modes in a plane-wave approximation. One of the motivations for the present 
paper is to provide a framework for computing the form of the magnetically modified global r- 
modes. The formalism and calculations presented in this paper are only for global oscillation 
modes of an infinitely conducting fluid. We have not included features such as a crust or 
a magnetosphere which would exist in the case of neutron stars. It should be possible to 
extend the present formalism to include interactions with crustal oscillation modes and a 
magnetosphere. 

The structure of this paper is as follows. In section 2 we review the equations describing 
the modes of rotating magnetic fluid stars. The form of the equilibrium magnetic fleld is 
discussed in section 3. In section 4 we show how the frequencies of the r-modes are modifled 
by the presence of a dipole magnetic field. We examine two model interior fields, the uniform 
magnetic field and the force-free magnetic field. In the case of a uniform magnetic field the 
frequency in the rotating frame of an r-mode with quantum number i is 



i + T 10^ ' 



sin^ a 

2 



(1) 



where Q is the angular velocity of the star, Ai is the energy stored in the magnetic field, T is 
the star's rotational kinetic energy and a is the angle between the magnetic field's axis and 
the spin axis. In the case of the force-free magnetic field, numerical results are presented. 
Finally, we conclude in section 5 with remarks about further problems which can be attacked 
with the present formalism. 



2. Normal Modes of Rotating Stars with Magnetic Fields 

In this paper we study the magneto-hydrodynamic (MHD) perturbations of a uniformly 
rotating star endowed with a magnetic field. As we do not assume that the symmetry axis 
of the magnetic field is aligned with the star's rotation axis, magnetic multipole radiation 
is emitted, and no equilibrium solution exists. However, we only consider timescales short 
compared to the characteristic stellar spin-down time due to the nonaxisymmetric magnetic 
field. In this sense an equilibrium solution exists. We work within the ideal MHD framework 



- 5 - 



(Jackson 1975) where the star's fluid is assumed to be electricaUy neutral and to have infinite 
conductivity. With the assumption of infinite conductivity, the magnetic field lines are frozen 
into the fluid, and the magnetic field pattern rotates at the same rate as the star. 

The equilibrium star is described by the density p, fluid pressure p, gravitational po- 
tential 0, velocity v and magnetic fleld B. Eulerian perturbations of these quantities are 
denoted with the symbol S. For inflnite conductivity, the Eulerian perturbations of the 
magnetic field and current density are 

= V X (5v X B) (2a) 

5E = --5vxB (2b) 
c 

-SJ = J- X (5B - -^Se] , (2c) 

c 47r V cdt J ' ^ ^ 

where = dt^ and the perturbed density, pressure and gravitational field arc given in 
standard works on perturbation theory such as Unno et al. (1989). The fluid displacement 
vector ^ is a solution of Euler's equation, 

'^ + 2Ck + yV^^ -F"^"^ , (3) 
P 

where the operators C and W are deflned by 

C^ = nx$, W^ = -V5p-^6p + V5(P, (4) 

P P 

and the magnetic force is 

pmag = 1 ((^J x B + J X 5B) . (5) 

The possible perturbations of a magnetized rotating star are found by specifying the equi- 
librium configuration of the star and then solving equation (3). This is a nontrivial problem, 
and the rest of this paper will be concerned with a method for finding normal mode solutions 
of Euler's equation. 

In this paper we adopt the method introduced by Schenk et al. (2002) where the per- 
turbations arc expanded in the spatial cigcnfunctions of Euler's equation with zero magnetic 
field. The normal modes of rotating stars with no magnetic fields (or any other external 
forces) are solutions of Euler's equation (3) with F"*"^ = of the form e~^^^*^A{x) where 
the spatial eigenfunction ^Ai^) satisfies the equation 

-u;\^A - '2tuJAC^A + >V£a = . (6) 
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Upper-case Latin subscripts are used to label solutions and correspond to the unique set 
of quantum numbers which describe the solution. We will make use of the inner product 
defined by 

< ^A, >= / d'xp{x)eA ■ C^B, (7) 

for any operator O. An important property of the modes of rotating stars discussed by 
Schenk et al. (2002) is that the modes are not orthogonal with respect to the usual inner 
product. Instead, the modes obey the relation 

Scarab = loa{uja+^b) < Ib, > +2'^^ < ^s, iC^A >= ^a^b < $a > +'^a < W|a > 

(8) 

which defines e^, the rotating frame energy of the mode (Friedman & Schutz 1978a). 

The lack of a simple orthogonality relation for the modes of rotating stars has the 
following consequence. If a general perturbation is expanded into a sum over modes of the 
form ^{x, t) = ciB(t)$B{x)j the expansion coefficients are not given by the usual formula 
aA{t) oc< ^A)C(^)^) > ^ would be true in the case of zero rotation. An alternative phase 
space expansion has been introduced by Schenk et al. (2002) to circumvent this complication. 
In the phase space approach, the perturbation and its first time derivative have expansions 
of the form 

^{x,t) = ^qB{t)^B{x), and ^{x,t) = -i^ujBqB{t)^B{x), (9) 

B B 

where the expansion coefficients are given by (Schenk et al. 2002) 

(t) = ^ {^A, u^A^ir, t) + t) + 2tC^{r, t)) . (10) 

The equations of motion for the expansion coefficients take on the simple form derived 
by Schenk et al. (2002) 

qA{t)+iuAqAt) = 1^ (^A, ^F-'>n^(r,^)]) = |f E (^A, -FsiUr^)]^ , (H) 

where the magnetic force (5) has an expansion of the form F"*"f = XIa ^^(^)-^>i[^a('')]- 

The term is the Lorcntz force created by the fluid motion We now deflne 

dimensionless magnetic coupling coefficients through 

KAB = ^< ^A, -Jb > ■ (12) 
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The integral < ^F^ > is the work done by the force when the fluid is displaced 
through the distance The coefficients kab can then be thought of as the ratio of work 
done by the perturbed Lorentz force to the total magnetic energy stored in the equilibrium 
star. The diagonal entries of the matrix k, correspond to the work done against the fluid 
motion which generates the Lorentz force F^. Since the magnetic forces will tend to 
oppose the motion of the fluid, we expect that the diagonal entries in the matrix of coupling 
coefiicients will be negative. 

After integrating by parts, the coupling coefficients take on the form 

-L_ J (^,« (,) - (,) - fn%ir)^ rHr (13) 

= y" V X (C^ X B) • V X (^^ X B) dQ, (13a) 

= y"(£:; X (V X B)) • V X X B) dn (13b) 

= J{CaX^)- (€d X B) dn. (13c) 

The diagonal elements /t^]^ and positive definite. In most cases ljR <^ c so that the 

term proportional to k^^, can be neglected. This last term is large typically only for neutron 

(3) 

stars. However, since we are using Newtonian gravity, the errors incurred in neglecting 
are similar to the errors coming from the neglect of general relativity. Although most authors 
neglect the term k^jdi '^^ have chosen to keep it explicitly in our equations for the sake of 
completeness. 

Normal mode solutions to the MHD perturbation equations can be found by making 
the ansatz 

QAit) = CAex.p{-iat) (14) 

where the ca are constants and a is the mode frequency in the rotating frame. Substituting 
the ansatz (14) into the equations of motion, we reduce the equations to the form 



M 

D 



oja\ca- — ^ kadCd I = crcA, (15) 



which is an eigenvalue problem for the frequency a and eigenvector components ca- In order 
to solve the eigenvalue problem, we must truncate the system to a finite number of modes. 
Suppose, for instance, that we chose to truncate the basis to a set of N eigenmodes of the 
non-magnetic system. Solving the system of equations (15) will result in N eigenvectors and 
eigenfrequencies of the magnetic system. The dimension of the vector space can then be 
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increased by one by adding one more non- magnetic basis vector and solving for the N+1 
magnetic modes. The difference between the original N frequencies can then be calculated 
and the process iterated until the frequency differences between successive iterations is suit- 
ably small. This procedure would be difficult to implement if it happened that all of the 
modes were coupled to each other. However, in the case of dipole magnetic fields, it will 
be shown in section 4 that the angular momentum selection rules restrict the possible mode 
couplings so that each mode couples to only a small number of other modes. As a result, 
most of the kad coupling coefficients vanish. This result will continue to hold for magnetic 
fields of higher mulitpoles. 

The constants appearing in the eigenvalue equation (15) are the rotating frame 
energies of the modes of a nonmagnetic star. These energies depend on the normalization 
chosen for the modes. Since the problem of finding modes of a magnetized star is linear, 
the final solution corresponding to an eigenfrequency a and eigenvector components ca is 
independent of the normalizations chosen for the spatial eigenf unctions of the nonmagnetized 
star. This allows us to choose any convenient normalization scheme which simplifies the 
form of the eigenvalue equation. A natural choice of normalization results by noting that the 
rotating frame energy is proportional to the square of the nonmagnetized mode frequency uja- 
The mode frequencies of rotating nonmagetized stars fall into two classes, frequencies which 
in the nonrotating limit are finite (such as f-modes, p-modes and g-modes) and frequencies 
which vanish in the nonrotating limit (inertial modes, such as the r-modes). A simple scheme 
is to choose the normalization of all modes with finite frequency in the nonrotating limit so 
that their energy is = where W is the gravitational potential energy of a constant 
density star with the same mass and radius as the star under study. Similarly, a suitable 
normalization is to choose the energy of the inertial modes to be = T where T is the 
rotational kinetic energy of the equivalent constant density star. This choice of normalization 
suggests the following approximation scheme for weak magnetic fields. When Ai/T <^ 1 and 
A^/ 1 I -C 1 the off-diagonal entries in the matrix 



will be much smaller than the diagonal entries. When Nad can be approximated as diago- 
nal, Nad — Sad ^1 — 2^^aa^ eigenvalue problem is greatly simplified, yielding spatial 
eigenfunctions identical to the eigenfunctions of the nonmagnetic star. In the case of modes 
which vanish in the non-rotating limit, such as the r-modes, only one mode exists for each 
spatial eigenf unction. In this case, when matrix Nad can be approximated by a diagonal 
matrix, the equation for the MHD frequency is simply given by 



Nad = Sad - - — k-ad 



(16) 




(17) 
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If the diagonal terms in the couphng matrix k are purely real, as in the case of a force-free 
magnetic field, the resulting frequency a will also be real. In the case of r-modes of stars 
with force-free magnetic fields, the term i^^'^j^{r) vanishes and the diagonal entries kaa are all 
negative definite. As a result, the magnetic field increases the absolute value of the r-mode 
frequencies. 

Implicit in this method is the assumption that the spatial cigcnf unctions of a rotating 
star form a complete basis. However, it has been shown by Dyson & Schutz (1979) that 
the normal mode of a rotating star arc complete only in a weak sense. Dyson & Schutz 
(1979) have shown that when the frequencies of the normal modes arc all real, there are 
no exponentially growing solutions of the initial-value problem. In the case of purely real 
frequencies it is then possible to express perturbations in terms of normal modes. Dyson 
& Schutz (1979) have stressed that their result allows the approximation of an infinite- 
dimensional problem by a finite-dimensional problem. These results on completenesss are 
only valid when the eigenf unctions are square intcgrable. In some special cases, such as 
an incompressible perfect fluid confined to a spherical shell, it has been shown (Rieutord 
et al. 2001) that there exists perturbations which are not square integrable and that only 
the r-modes are regular solutions. For the case of a spherical shell, an expansion in normal 
modes will not represent all perturbations. 

Clearly, this method relies heavily on the computation of the magnetic coupling coef- 
ficients K,A£). This is not a simple task. In the following sections we discuss selection rules 
and strategies for computing the coupling coefficients. 



The calculation of the kad coupling coefficients is nontrivial and requires a systematic 
strategy. In our computation we make use of spin-weighted spherical harmonics and the 
formalism introduced by Schenk et al. (2002). 

To begin with, we introduce a complex basis {cq, e+, e_} related to the usual orthonor- 
mal basis {e^, e^, e:} by 



2.1. Spin-weighted Mode Decomposition 




(18) 



These unit vectors have the following scalar products 



Co • Co = e+ ■ e 



1, 



eo • e± = e+ • e+ 



= e_ • e_ = 



0. 



(19) 



Any vector ^ has an expansion of the form 



(20a) 
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Co^Cr, e± = ^-e± = -^(ee±^e^), (20b) 
where Ce, C<t> the components of the vector in the usual orthonormal basis. 
The spherical harmonics with spin-weight s are defined by (Campbell 1971) 



Xem{9, 0) = ^^^dU^{e)e^"^f (21) 

where the dlj^{9) are the matrix representations for rotations through an angle 9 discussed 
in detail by Edmonds (1974). When the spin- weight s = 0, the spin- weighted spherical har- 
monics reduce to the regular spherical harmonics. Prom the definitions of the spin spherical 
harmonics, it follows that (Schenk et al. 2002) 

e±-VYim = T^VW+^iiYem (22a) 
v2r 

e±-LYem = -^^£{£+1) ±1^,^, (22b) 

where L = —ir x V. 

The spatial eigenmodes of a rotating, unmagnetized star have azimuthal angular depen- 
dence exp(imyi0), but in general do not have a definite quantum number £. The modes can 
be written in the general form (see, for example Lockitch & Friedman (1999)) 



e=\mA\ 



(23) 



r r 
Making use of equations (22a) and (22b), the general mode expansion (23) is 

= E {C^'W^mAeo + fr^{r)+,Y,^,e. + f_"'- {r) _,Y,^,e+) (24) 
i' 

where the functions /f"*^ are given by 



ciruA ( 



f^^^ir)^:!^ and fr\r) ^ ^\r-^^-{V,mA^U,^,) . (25) 



3. The Equilibrium Magnetic Field 

A static magnetic field which is symmetric about an axis can be written in the general 

form 

B(r, 9\ 00 = E {Mr)yto{9\ 0')^. - r/32,^(r) VY,o(^', 0') - i^^,^{t)VY,o{9\ 0')) , 

(26) 



- 11 - 



where 9' and 0' are angular coordinates measured with respect to the magnetic field's axis of 
symmetry and V and L' are angular operators with respect to the magnetic field's symmetry 
axis. The expansion of the general magnetic field in the spin- weighted basis is 

S(r, e\ <P') = B',{r, 9', 0')eo + B'^{r, 9', cj>')e_ + B'_ir, 9', 0')e+, (27) 

where each component has the expansion in spin-weighted spherical harmonics 

B's = Y.b'lir)^eoie', (28) 
e 

and the set of functions b'l{r) are given by 

'"'-^A' ^l^i/^^/^Wv^'A.). (29) 

In general, the axis of symmetry of the magnetic field will not coincide with the star's 
spin axis. Suppose that the magnetic field is tilted by an angle a from the star's spin axis. 
If the angles 9 and arc measured with respect to the spin axis, the spin-weighted spherical 
harmonics in the two coordinate systems arc related by 

rek{9', 00 = "^'ruM ^UO, 0). (30) 

m 

As a result, the magnetic field can also be written in terms of the spin-axis coordinate system 

5(r,^,0) = Boir,9,(t>)eo + 5+(r,^^,0)e_ + 5_(r,^,0)e+ (31a) 
B, = J2b'rir)sYem{0:<P) (31b) 

£,m 

b'rir)=b'l{r)dUo^). (31c) 

In order for the magnetic field to be real, the coefficients must obey the relations 

bTir) = b'rir) (32a) 
6j™*(r) = -b'f'ir) (32b) 

where * denotes complex conjugation. 

The exterior magnetic field of most stars is approximately that of a dipole, which has 
the form B(r', 0') = -BpR^y/An/3 V{Yio{0' ,(/)')/ r'"^) in the coordinate system aligned 
with the magnetic field's symmetry axis. In this case the exterior magnetic is described by 
coefficients Pi{r)/2 — P2{r) — BpB?/r^ and I3z{f) — 0. (We will now drop the subscript I 
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for the P coefficients since we will only be considering the dipole case £ — 1.) In the case of 
a dipole field, the only coefficients of the spin matrix d}{a) of interest are 

^io(<^) — — ^sina, (ioo(<^) ~ cos a, d\Q{a) — 7= sin a. (33) 

In the coordinate system aligned with the star's spin axis, the exterior magnetic field is then 
described by the functions 

-bl"\r) = hl-{r) = -6i-(r) = \j^dUa)B,-. (34) 

In order to compute the eigenmodes of a star with an dipole magnetic field exterior to 
the star, the magnetic field inside of the star must be specified. In this paper we will only 
consider interior fields which are dipole (i.e. i — 1). General formulae for the perturbations 
of a star with any interior dipole field will be presented in the following section. As examples 
we will consider two simple interior magnetic fields which can be matched to the dipole 
exterior (in the sense that the normal component of the magnetic field at the star's surface 
is continuous) . Both of these examples are force-free fields which obey 

V X B = /xB, (35) 

where /i is a constant with units of inverse length. Force-free magnetic fields have the simple 
feature that they do not distort the star's equilibrium fiuid configuration. 

The first example is the simplest possible interior field, the uniform magnetic field, 
corresponding to // = in equation (35). The coefficients appearing in equation (26) are 

/3i = -/52 = 5,„ /53 = 0, (36) 

where 5j„ is the magnetic field at the centre of the star. Since the radial component of the 
magnetic field must be continuous at the surface of the star. Bin — '^Bp. The magnetic field 
coefficients in the coordinate system aligned with the spin axis are 

b'o-'ir) = bh-'ir) = -b'rir) = ^B.^dU^). (37) 

The second example corresponds to the force-free condition (35) with /j, nonzero dis- 
cussed by Ferraro & Plumpton (1966). The magnetic field coefficients are 

P,{r) = 3Bj-^, /32(r) = -^B,„i(^|; + i^ji(H, /Ss = -^i^^nJi(H, (38) 
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where ji{nr) is a spherical Bessel function. Continuity of the normal component of the 
magnetic field at the surface of star requires that the constant be given by the solution of 
the equation 

hiliR) = I^IJ^R, (39) 

which can be solved numerically once the ratio Bp/ Bin is specified. The solution to equation 
(39) is shown in Figure 1, where the dimensionless product is plotted versus the ratio 
Bp/ Bin- The force-free magnetic field then corresponds to family of fields parametrised by 
the ratio of magnetic field at the pole to the magnetic field at the centre of the star. This 
ratio of magnetic fields can have values in the range < Bp/ Bin < 0.5, as shown in Figure 
1. 



4. Magnetically Modified r- modes 

We now turn to the calculation of the modification of the r- modes due to dipole magnetic 
field. In the slow-rotation approximation, the r-modes of non-magnetic incompressible stars 
have the form 

Wt^ruA = 0, Vt^m^ = (40b) 
UiAm. = K\mAR y j (40c) 

e^^m^ = T = \MR^n\ (40d) 
o 

where the functions We^mA^^iAmA U^^m^ are defined in the expansion (23). Each r- 
mode solution is identified by a distinct quantum number itla- At lowest order in angular 
velocity, the only value of I a allowed is I a = I^aI- The spin- weighted decomposition of the 
Aih r-mode is 

= foir) ^tAruA^O + f+{r) ^iAruA^- + f-{r) ^iAruA^^. (41) 

where the functions f^{r) are given by 



/„V) = 0, /^ = /_- = fl(«, + l)^^P^±±^(L)'\ (42) 



In order to compute the r-mode frequencies in the presence of a magnetic field, we must 
calculate the values of the magnetic coupling coefficients appearing in equation (13). The first 
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coupling coefficient K^Aoir) depends on the perturbed magnetic field. The equations for the 
perturbed magnetic field were derived in the appendix for general magnetic fields and general 
fiuid perturbations. In this section we will present the form of the perturbed magnetic field 
for r-modes in the presence of a dipole magnetic field. In general, the perturbed magnetic 
field has the expansion 

SBA{r, e,4>) = Y. i^^Aoi^) o^veo + SB'^^^t) _y,;e_ + SB'^^_{r) ^,;e+) . (43) 

In this equation, the summation is over all A obeying the triangle inequality i?A — 1 < A < 
iA + i and all /i obeying fj, — —Ia — m. In the case of an r-mode and a dipole magnetic field, 
the "0" component appearing in equation (43) is given by 



m=-l ^ 



^ ')(y'\ ]) -b'rm (1 + (-1)^+^-) , (44) 

II ruA m/\0— 11/r ^ ^ 



which is only nonzero when = —{mA + m) and \ — Ia and the constant C(^a, 1, A) is given 

by 



(:;(£^,£,A) = (-i)^-+^+^ 

The ± components appearing in (43) are 

= (-1)^+^A+1^5A,M* 

= C{iA. 1, A) [ 

m=-l ^ 



(2£^ + l)(2£+l)(2A + l) 
47r 



A I, 



' ) -fiir) X 
/X TUa m J r 



{{iA - l)&r (0 - bTir)) 



A iA 1 \ ,/ ^a(^a + 1) , 
-1 1 / V 2 + ^ ^ 



{iA-mA+2) ^,,^^^^ 



(45) 



(46a) 
(46b) 



A iA 1 
-1 1 



A £a 1 
-1 2 -1 



In equation (46b) A takes on values iA,iA ± 1- Evaluating the terms appearing in square 
brackets in equation (46b), the individual terms for the allowed values of A are 



(-1) 



3(2£^ + l) f^jr) 
AnlAiiA + 1) r 



X 



E 

m=— 1 



tA - iwr - {bT + bT*) - i^UiA + i)(&i' 



l.m 



7 l,m* \ 

b+ ) 



-{iA + m) £a m 
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(-1) 



l3 {£A-mA + l)ft{r) 



X 



E 



m=— 1 



l,m I 7,1, m*> 



+1 



m=— 1 



Ia-I iA I 
-{iA + m) iA m 



iA + l iA I 

-{iA + m) iA m 



The first coupling coefficient k^]) is given by equation (A13). The first summation 
appearing in equation (A13) reduces to 



J2SBto'5B'^'o-^SADC'{iAAjA) i: 



m=—l 



iA iA 1 

H ruA m 



iA iA 1 

-1 1 



(48) 

Making use of the known values of the Wigner 3-j symbols (Edmonds 1974), equation (48) 
reduces to 



iA + 



-Sad (^)) (^^A cos' a + ^ sin' . (49) 



Due to the symmetry property (46a), the final two summations appearing in (A13) are 
^ (SB*^'XSB'^% + 5BXi'SB'^^_) = Yl ^Kt^B^n'-^ + {-iY^^'''SB'x^_^SB*^';^. (50) 



Since the triangle inequahties ^a — 1<A<£a + 1 and ^d — 1<A<^£) + 1 must both be 
satisfied, only modes satisfying i^ — iA,iA + i,iA + '^ have non-zero magnetic coupling. As 
a result the only nonzero entries in the coupling coefficient matric kad will be those coupling 
the modes satisfying the triangle inequality. 

The general equation for the diagonal elements of the coupling coefficient matrix K^]i(r), 
for any interior dipole magnetic field is 



4'i(r) 



2 ('Mi l X 



r—T-riPif i^A cos' a + J sin^ a 
'A + J- V ^ 



(51) 



iA{i. 



{{{iA - m - 2/^2)' + {Psf) (^iA cos' a + i sin2 
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2£a(2£a + 1) 
Similarly, the function k^]^ has the value 



m - m - {^A + 2)hf sin^ a 



{2i\ + 3iA + 2) 
2(2£a + 1) 



- ( 



U{r) 



(£a + 1)(2£a + 1) 



X 



d 



2/3i ( ^(r/?2) + A ) - 2 {PsY ) ( £a cos^ a+- sin' a 



(52) 



The frequency correction can now be computed for any interior dipole magnetic field, once the 
functions l3i{r),P2{r) and Psi^r) which describe the equilibrium magnetic field (see equation 
(26)) have been determined. 



4.1. Frequency Correction for a Uniform Magnetic Field 

The simplest example of an interior magnetic field is the uniform magnetic field which 
has f3i{r) = —l32{r) = Bin and ^^{r) = 0. Since the magnetic field is uniform, it satisfies 

(2) (3) 

the equation V x B = 0. As a result, the term k)^^ vanishes. The term does not 
vanish, but its contribution to the complete coupling coefficient kad is a factor of Vt^R^ /c^ 
smaller than the contribution from h^adi neglect it. The magnetic coupling coefficients 

between different r-modes can now be computed. The only nonvanishing coefficient is the 
self-coupling term kaa given by equation (51). A straight-forward calculation yields 



«AA = -^(^A+l)(2^A + 3) 

o 



l + ^(^A(^A + l)-3)sin2Q; 



(53) 



Since the matrix k is diagonal for the case of coupled r-modes, the frequencies cr^ of the 
magnetically modified r-modes are simply 

= '^A (^1 + ^I'tAAl^ ■ (54) 

The final equations (53) and (54) for the frequency correction for the r-modes is sur- 
prisingly simple. However, this simplicity is a result of the simplicity of the background 
magnetic field. Since we assume that the field is homogeneous, terms involving derivatives 
of the magnetic field vanish, with the net result that the matrix had is diagonal. 
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In this simple example, we have not considered couplings between r-modes and other 
families of modes, such as the f-modes or inertial modes. By neglecting these couplings, we 
are neglecting off-diagonal terms in the matrix Kadi which is only valid in the approxima- 
tion that the energy in the magnetic field is small compared to the star's rotational and 
gravitational potential energies. For this reason equation (54) is the frequency correction 
to the r-modes only in the weak-field approximation. We will compute the couplings to 
other modes elsewhere, which will allow the calculation of the modes of rotating stars with 
arbitrarily large magnetic field. 



4.2. Frequency Corrections for a force-free Magnetic Field 

The second model for the interior magnetic field which we have examined is the force- free 
magnetic field, which is described by the functions /3i(r), /92(r) and jSzij-) given by equation 
(36). All three functions arc proportional to which is the magnetic field strength at 
the centre of the star, which is not a known quantity. Typically, only the exterior dipole 
magnetic field at the pole of the star. Bp, can be measured, and the interior magnetic field 
is unknown. When the interior field is modeled by a force-free field, the ratio of Bp/ Bin is 
a free parameter. The smallest ratio of interior to exterior magnetic fields allowed in the 
force-free magnetic field model is Bp/ B^^ = 0.5, which corresponds to the uniform magnetic 
field (see Figure 1). 

Due to the toroidal symmetry of the fiuid displacement vector, k/^j^ vanishes for the 
force-free magnetic field so that only needs to be computed. This computation can be 
easily done numerically and the results are shown in Table 2 and Figure 2. The frequency 
correction depends on the strength of the internal magnetic field and the angle between 
the spin axis and magnetic field axis. In order to show the dependence on both of these 
parameters, we have written the correction term in the form kaa — claa + ^aa sin^ a, where 
aAA and hAA depend only on the ratio of external to internal magnetic field Bp/ Bin- In Table 
2 we show the values of the magnetic frequency correction terms for various values of Ia and 
for the specific value of magnetic field ratio Bp/B^^ = 0.1. It can be seen that the couphng 
coefficients grow with increasing values of the spherical harmonic index Ia- This property 
agrees with what is seen in the case of the uniform magnetic field. 

In Figure 2 we show the values of the frequency correction for the case Ia — tua — 2 
for the complete range of allowed interior magnetic field strengths. Figure 2 shows that the 
magnitude of the magnetic coupling coefficients decreases as the interior magnetic field is 
increased, which seems counter-intuitive. However, in order to find the corrected frequency 
in equation (54), the coupling coefficient is multiplied by M., the energy of the internal 
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magnetic field. Now M. oc Bf^, and the internal magnetic field is not known: all that 
we can measure is the external dipole magnetic field strength Bp. Once Bp is known, the 
fraction Bp/ Bin is a free parameter of the theory. Therefore, the true scaling of the frequency 
correction due to magnetic fields is {Bin/ Bp^KAA- Once this is taken into account, we see 
that the magnitude of the frequency correction increases with increasing interior magnetic 
field strength, as would be expected. 

5. Conclusions 

In this paper we have presented a formalism for computing the oscillation modes of 
magnetic rotating stars. This method is based on the phase-space mode decomposition 
introduced by Schenk et al. (2002) and makes use of spin-weighted spherical harmonics. 
As a simple example, we have computed the magnetic corrections to the r-modes of an 
incompressible star with a uniform magnetic field in the weak-magnetic field limit. In this 
limit, the spatial form of the r-mode velocity field is unchanged by the presence of a magnetic 
field. In the case of a uniform interior magnetic field matched to an exterior dipole field, the 
resulting formula for the frequency has the simple form presented in equations (53) and (54). 
Qualitatively similar results were found for the force-free magnetic field. The presence of a 
magnetic field increases the absolute value of the rotating frame frequency, so that the mode 
will be counter-rotating faster than if the star had no magnetic field. Adding a magnetic 
field to the star increases the restoring force acting on the fiuid and increases the oscillation 
frequency, just as increasing the spring constant of a spring causes a spring to undergo faster 
oscillations. Mathematically, the increase in frequency results from the fact that the diagonal 
entries in the matrix k are negative. Since these diagonal entries correspond to the work 
done by the magnetic field on the fluid, the negative sign denotes the fact that the magnetic 
field opposes the fluid motion. Hence, the increase in frequency which we have found agrees 
with basic physical intuition. 

Although we have only presented detailed results for two types of internal dipole field, 
the equations presented in section 4 can be used to determine the frequency correction for 
any other internal dipole field. Once the functions /9i(r) — (^^^{r) describing the equilibrium 
magnetic field have been found, it is only necessary to substitute them into equations (51) 
and (52) and perform the integration over volume in equation (13). 

As long as the magnetic field energy is small compared to the star's rotational energy, 
the CFS gravitational-radiation-driven instability criterion Friedman & Schutz (1978b) will 
be unchanged. As a result, the magnetic field has a slight stabilizing effect on the r-modes, 
since it is now harder for the mode to be dragged forward by the star's rotation. This 
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agrees with the intuitive notion that magnetic fields tend to oppose fiuid motion. However, 
in order for the mode to be stable to the CFS mechanism, it would be necessary for the 
magnetic field's energy to be of the same order of magnitude as the rotational energy, which 
corresponds to a regime where our present results do no hold. Such a large field would 
presumably alter the stability criterion (Friedman & Schutz 1978b) as well. 

Although our calculations for the r-modes of magnetic stars is only valid in the weak- 
field limit, we note that a young neutron star with spin period and magnetic field similar to 
the Crab's will have M./T of the order 10~^ well within the regime of the weak-field limit. 
Since we do not include a crust, the applicability of our equations is mainly to young neutron 
stars with T > lO^^K or other stars without a solid layer. 

Our calculations are purely linear, so the present work does not shed any light on the 
nonlinear evolution of an unstable mode of a magnetic star which has been investigated by 
RezzoUa et al. (2000). However, the frequency change due to the magnetic field which we 
have calculated will alter the kinematic drift found by RezzoUa et al. (2000), although we 
doubt that these changes would alter their qualitative results. It would be interesting to see 
how the magnetically modified r-modes presented in this paper affect the evolution of spin 
and magnetic field in young neutron stars. 

This research was supported by the Natural Sciences and Engineering Research Council 
of Canada. We would hke to thank Roy Maartens, Draza Markovic and Luciano RezzoUa 
for useful discussions and comments about this paper. 
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A. Derivation of the Magnetic Coupling Coefficients for General Magnetic 

Fields 

In this appendix we give details about the calculation of the magnetic coupling coef- 
ficients which appear in equation (13). The first coupling coefficient /t^]) depends on the 
perturbed magnetic field, defined by = V x [C,a x B), where B is the equilibrium mag- 
netic field and is the fiuid displacement vector. This definition of the perturbed magnetic 
field can be rewritten in the form 

dB^^ = {aB'-eAB''h, (Ai) 

where the subscripted semicolon denotes covariant differentiation, and lower case Latin sub- 
scripts and superscripts arc vector indices. The components of the perturbed magnetic field 
in the basis (eo,e+,e_) can be found by substituting the expansions (20a) and (31a) into 
equation (Al). The radial component of the perturbed magnetic field is given by 

55° = SBA-eo^eA;,B' + eAB\-eA;,B'-eAB% 

= a ■,^B^ + ^A .-B- + + - ^A - a - ei^v - a5fA2) 

This preprint was prepared with the AAS IM^iX macros v5.0. 
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Note that although it is true that B''.^^ = 0, we gain an equation with a more symmetric form 
by not exphcitly setting the divergence of the magnetic field to zero. Since each component 
of the equilibrium magnetic field and the fluid displacement vector is proportional to a spin- 
spherical harmonic, the perturbed magnetic field components can be written in the compact 
form 

6BAAr,0,<i^) = J2Yl E CLsir) s.Ye,^, sYem (A3) 

£a t,m sa+s=s 

where s takes on values of and ±1, and the radial functions C^i^,s(^) will be determined 
shortly. In equation (A3), the integers iAjiriAjSA are the fluid perturbation's quantum 
numbers, i, m, s are the equilibrium magnetic field's quantum numbers and the symbol A 
represents the set of quantum numbers {£, m) describing the magnetic field. 

A detailed formalism for taking covariant derivatives of vectors using spherical coor- 
dinates has been given in section VI of Schenk et al. (2002). Their equations (6.15) and 
(6.16) summarise the rules for taking covariant derivatives in the spin- weighted basis. (The 
only difference in notation between the present paper and that of Schenk ct al. (2002) is the 
symbols used for the basis vectors. Schenk ct al. (2002) use the basis {I, m, rn} which are 
related to our basis by I = eg, m = e+ and m = e_.) When these covariant differentiation 
rules are used, we find that the s = terms of ("^^ appearing in (A3) are 



aUr) = a/^^^(&^-&^)-a/^^^%^^(/,"-/-") (A4a) 



r 




Similar manipulations give us the ± components of the perturbed magnetic field. In 
particular, the -|- component of SB is 

SB+ = SB- = SB-e+= C.tB' + - C^^" - ^B', 

-e,+B- - e,oB- - eB-,^ - eB-,o . (as) 

Substitution of the covariant differentiation rules into (A5) and comparing with equation 
(A3) results in 



-23- 



In the case of a purely dipole magnetic field, £ — 1 and the term Ci^:!^i,2(^) vanishes. Similarly, 
the s = — 1 components of ( are 



et.(.) ^ /^ZIfZ±5^ (Arc) 



A fA 



In equation (A3) , the components of perturbed magnetic field were written as products 
of spin-spherical harmonics. This equation can be simplified so that each term is proportional 
to only one spin-spherical harmonic, 

5B^,,(r, 9,<P) = J2 SBZir) X^,, (A8) 

by using the combination formula 

^^ ^^ V- /(2£ + 1)(2A; + 1)(2A + 1) f £ k X \ f £ k X\ ,^ 



(A9) 
< A < 



where the summation is over all values of A satisfying the triangle inequality \l — k\ _ 

( £ k A \ 

£ + k, fj, satisfying fj, — —{m + n), a satisfying a — —{s + t) and I ) 

y —s —t —a J 

Wigner 3-j symbol (Edmonds 1974). The radial expansion functions SB^^{r) are given by 
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/ 



and the constant C{£a,£, A) is defined by 

cv., e, X) = ^_,y...., fWlTmi+im+i)^ ^^^^^ 

V 4:71 

Since the spin-spherical harmonics obey the orthogonahty relation 

'47r 

the first magnetic coupling coefficient is just the summation 

>^%ir) = [SB*^. 6Bndn = {sBf/6B'^^, + 6BfX6B'^^_^ + 6Bf;_^6B'^^_) . (A13) 

Similar expressions for the coupling coefficients K^])(r) and K^])(r) can also be written. 
In order to calculate K^])(r), the vector ^ x (V x B) needs to be decomposed into a sum 
over the spin- weighted spherical harmonics. The decomposition is of the form 

= X (V X B) = 5] (xlfo(r) o^^eo + x^A^) -Y^.e. + xt-i^) ^a>+) , (A14) 
\,n 

where the coefficients are 

The only non-zero functions x'^^(r) are given by the expressions 

= (^0 + , Xo^ti(r) = - (r6^),.) (Al6a) 

xio,i(r) = - J/o^(6^ + K),.) , xi;to(r) = - + h^-) (Al6b) 

X^to,-i(0 = Vcf (fc^ - (r6^),.) , W = + &^)- (Al6c) 
The second coupling coefficient K]^jy is now given by 

^%{r) =Jxa- SBndn = {xto'^B'^^o + Xa'X^B'^'-, + X^l^^i^^'r-) • (A17) 

Similarly, we define the vector 

1a = ^axB = J2 <{fAo(^) o^A;eo + 7lr+(r) -r,;e_ + 7Af-(0 ^A;e+) , (A18) 
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where the expansion coefficients are 



Tfie non-zero coefficients 7'^^(r) fiave tfie values 



7o^^,i(r) = ^^^h% , ^^l_,{r) = -^fX- (A20a) 

7^^o(r) = , ^tu^) = -^ux (A20b) 

^-U,o{^) = -^/-^o , 7^to,-i(^) = -ifobt (A20c) 



Tlie tliird coupling coefficient is 
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Table 1. Typical ratios of magnetic field energy to gravitational and rotational energies 
for different types of stars, assuming uniform density and magnetic field. The energy in the 
magnetic field is A4 = B^^B? /Q, the gravitational potential energy is \W\ = 3GM'^/5R and 
the rotational energy is T = Mi?^(27r/P)^/5. White dwarf data is from Wickramasinghe & 

Ferrario (2000), Ro AP star data is from Kurtz (1990), and the data for SGR 0526-66 is 

from Duncan (1998). 



Name 


lype 


B 


Period 


Mass 


Radius 


M/\W\ 


M/T 






(G) 


(s) 


(Me) 


(cm) 






Crab 


NS 


l.Oe+12 


3.3e-02 


1.4 


l.Oe+06 


5.4e-13 


8.3e-09 


AP stars 


RoAp 


l.Oe+03 


l.Oe+06 


2.0 


1.4e+ll 


l.Oe-10 


7.4e-07 


AM Her 


WD 


1.3e+07 


l.le+04 


0.7 


l.Oe+09 


3.6e-10 


3.2e-04 


PG 1015+14 


WD 


1.6e+08 


5.9e+03 


0.7 


l.Oe+09 


5.5e-08 


1.4e-02 


PG 1031+234 


WD 


l.Oe+09 


1.2e+04 


0.7 


l.Oe+09 


2.1e-06 


2.3e+00 


SGR 0526-66 


NS 


l.Oe+15 


S.Oe+OO 


1.4 


l.Oe+06 


5.4e-07 


4.9e+02 



Table 2. The numerical values of kaa for force-free magnetic field. The data are 
calculated for fj, — 3.51835 and -f^ = 0.1. The coupling coefficients are written in the form 
Kaa — ciAA + ^AA sin^ a where a is the angle between the star's spin axis and the magnetic 

field's symmetry axis. 



^A 


dAA 


hAA 


1 


-4.2182e-01 


2.1091e-01 


2 


-5.9812e-01 


2.5302e-02 


3 


-1.5734e+00 


7.8447e-03 


4 


-4.1159e+00 


2.9866e-01 


5 


-9.2792e+00 


1.0233e+00 


6 


-1.8367e+01 


2.2962e+00 


7 


-3.2933e+01 


4.2251e+00 


8 


-5.4780e+01 


6.9128e+00 


9 


-8.5963e+01 


1.0458e+01 


10 


-1.288e+02 


1.4955e+01 
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Fig. 1. — Plot of IJ,R versus Bp/ Bin, for the force-free magnetic field. In this plot R is the 
radius of the star, /i is the proportionality constant appearing in equation (35) and Bp/ B^n is 
the ratio of external to internal magnetic fields. This plot was made by numerically finding 
the solution of equation (39). 
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Fig. 2. — The magnetic frequency correction for an ^A = = 2 r-mode as a function of 
Bp/ Bin, for the force-free magnetic field. The frequency correction is written in the form 
K'AA — oiAA + ^AA sin^ «, and the constants 022 and 622 are plotted. The solid curve is the 
graph of the function 022 and the dot-dashed line is the graph of 622- 



